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Two flavors of dynamical quarks on anisotropic lattices 



m 
o 
o 

(N 



(N 

> 

(N 

o 

O 
m 
O . 

<D . 
43 ■ 



X 



T. Umeda*, 1 S. Aoki, 2 M. Fukugita, 3 K.-I. Ishikawa, 1 N. Ishizuka, 1 - 2 Y. Iwasaki, 1 ' 2 K. Kanaya, 2 
Y. Kuramashi, 4 V.I. Lesk, 1 Y. Namekawa, 2 M. Okawa, 5 Y. Taniguchi, 2 A. Ukawa, 1 ' 2 and T. Yoshie 1 ' 2 

(CP-PACS Collaboration) 

1 Center for Computational Physics, University of Tsukuba, Tsukuba, Ibaraki 305-8577, Japan 
Institute of Physics, University of Tsukuba, Tsukuba, Ibaraki 305-8571, Japan 
3 Institute for Cosmic Ray Research, University of Tokyo, Kashiwa, Chiba 277-8582, Japan 
4 High Energy Accelerator Research Organization (KEK), Tsukuba, Ibaraki 305-0801, Japan 
'Department of Physics, Hiroshima University, Higashi- Hiroshima, Hiroshima 739-8526, Japan 

(Dated: February 1, 2008) 

We report on our study of two- flavor full QCD on anisotropic lattices using 0(a)-improved Wilson 
quarks coupled with an RG-improved glue. The bare gauge and quark anisotropies corresponding to 
the renormalized anisotropy £ = a s /at = 2 are determined as functions of /3 and «, which covers the 
region of spatial lattice spacings a a ~ 0.28-0.16 fm and rnps/mv ~ 0.6-0.9. The calibrations of the 
bare anisotropies are performed with the Wilson loop and the meson dispersion relation at 4 lattice 
cutoffs and 5-6 quark masses. Using the calibration results we calculate the meson mass spectrum 
and the Sommer scale ro- We confirm that the values of ro calculated for the calibration using 
pseudo scalar and vector meson energy momentum dispersion relation coincide in the continuum 
limit within errors. This work serves to lay ground toward studies of heavy quark systems and 
thermodynamics of QCD including the extraction of the equation of state in the continuum limit 
using Wilson-type quark actions. 

PACS numbers: 12.38.Gc 



I. INTRODUCTION 

In spite of recent progress in computer technology and 
numerical algorithms, extraction of continuum properties 
from lattice QCD remains to be challenging when dy- 
namical quarks are included due to large computational 
demands. One of methods for alleviating the difficulty 
is to improve the lattice action for a faster approach to 
the continuum limit. This has enabled us to carry out 
the first systematic extrapolation to the chiral and con- 
tinuum limits for the light hadron spectrum 0, . 

Another method which is effective for several quanti- 
ties is to introduce a space-time anisotropy. In Ref. 
we have shown that using anisotropic lattices with a 
larger temporal cutoff is efficient for reducing lattice ar- 
tifacts in thermal QCD, and have carried out the first 
well-controlled continuum extrapolation of the equation 
of state in quenched QCD. In finite temperature QCD, 
anisotropic lattices have been employed in the quenched 
approximation also to study transport coefficients 
pole masses 0, glueballs Q, and spectral functions 
Q H , where anisotropy was introduced to obtain more 
data points for temporal correlation functions. At zero 
temperature, anisotropic lattices have been employed to 
study charmonium states 0, 0, 0, glueballs [l5| . 
heavy hybrids 0,0, and also the pion scattering length 

E 

In this paper, we calculate the anisotropy parameters 
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for an improved full QCD action to contribute toward a 
systematic study of QCD with heavy quarks and at fi- 
nite temperatures. The calculation of anisotropy param- 
eters is not a simple task in full QCD because a couple 
of bare parameters have to be simultaneously adjusted 
to achieve a consistent renormalized anisotropy in phys- 
ical observables for both quarks and gluons. This tun- 
ing of bare anisotropy parameters is called "calibration" 
[r7j . We study two-flavor full QCD with an RG-improved 
gauge action and a clover-improved Wilson quark action, 
extending the combination of improved actions adopted 
by the CP-PACS Collaboration to anisotropic lattices. 
Carrying out simulations at several values of bare pa- 
rameters, we perform the calibration to determine the 
bare anisotropy parameters for a given value of the renor- 
malized anisotropy £ = a s /at as functions of the gauge 
coupling and bare quark mass. We study the range of 
parameters corresponding to a s m 0.28-0.16 fm for the 
spatial lattice spacing and mps/mv ~ 0.6-0.9 for the 
ratio of pseudoscalar and vector meson masses. Based 
on our previous study of finite temperature QCD Q , we 
concentrate on the case £ = 2 in this paper. 

Different choices of observables for the calibration will 
lead to 0(a) differences in the calibration results. We 
study this issue by comparing the results from two dif- 
ferent observables — pseudoscalar and vector meson dis- 
persion relations. We anticipate that the results of differ- 
ent calibrations will be useful for checking the stability 
of continuum extrapolations. As a test of the idea, we 
also perform a continuum extrapolation of the Sommer 
scale at £ = 2, by interpolating our measurement results 
to the calibrated points. 

This paper is organized as follows: We define our 
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anisotropic lattice action in Sec.[nJ and discuss our choice 
of £ = 2 and simulation parameters in Sec. 11111 The cal- 
ibration procedure is described in Sec. IIVI Results of 
two calibrations for £ = 2, using pseudoscalar and vector 
meson dispersion relations, are summarized in Sec. W\ 
Finally, in Sect. l*VTl we interpolate the measurement re- 
sults to £ = 2 to study basic properties of our lattices. 
We also test how the difference in the calibration affects 
physical observables. Sec. IVIII is devoted to conclusions 
and discussions. An appendix is added to compare our 
calibration procedure with another method based on the 
ratio of screening and temporal masses. 



II. ANISOTROPIC LATTICE ACTION 

We study full QCD with two flavors of degenerate light 
quarks. On isotropic lattices, we have made a series of 
systematic studies adopting a clover-improved quark ac- 
tion coupled with an RG-improved glue Q, U 0, 0, 
I20L l2ll |22|. In these studies, the clover coefficient csw 
was set to the tadpole- improved value using the plaque- 
tte in one-loop perturbation theory for the mean-field. 
This choice was based on the observations that the one- 
loop plaquette reproduces the actual plaquette expecta- 
tion values within 8% for the range of parameters stud- 
ied, and that the resulting value of the clover coefficient 
agrees well with actual one-loop value of it. This action 
was shown to give both a good rotational symmetry of 
the heavy quark potential and a small scale violations in 
the light hadron spectra at moderate lattice spacings [23| . 
At zero temperature, these good properties enabled us to 
carry out first systematic chiral and continuum extrap- 
olations of light hadron spectra and light quark masses 
0, Q. At finite temperatures, this combination of ac- 
tions was shown to reproduce the expected 0(4) scal- 
ing around the two-flavor chiral transition point [181 [24j , 
and was adopted in the first systematic calculation of the 
equation of state in lattice QCD with Wilson-type quarks 
[20|. Here, we extend the study to anisotropic lattices. 



A. RG-improved gauge action on anisotropic 
lattice 



On isotropic lattices, the RG-improved gauge action by 
Iwasaki |25j consists of plaquettes and 1x2 rectangular 
loops. Extending it to anisotropic lattices, the general 
form of the action is given by 

Sg = & I — E 0=8^ (*) + c i(%M + %(*))} 

[ lG x,i>j 

+7gX!( c o^4W + c\R m (x) + c' 2 R 4k (x)) \, (1) 



where i, j, k are for spatial directions and 



Rnv{x) 



l--ReTr{U fl (x)U v (x + fi) 
xUt(x + uM(x)}, (2) 
1 - -ReTriU^U^x + fi)U u {x + 2/t) 
xUUx + fi + u)UUx + tk)Ut{x)}, (3) 



are plaquette and rectangular loop in the \irV plane, re- 
spectively. The improvement coefficients cf satisfy the 
normalization conditions, cg+8cf = 1 and Cq+4c\+ 4c| = 
1. The bare gauge coupling equals (3 = G/g 2 , and 70 rep- 
resents the bare anisotropy. We have three independent 
improvement parameters among cj l . 

In principle, these improvement coefficients may have 
non-trivial £ dependences depending on the improve- 
ment conditions on the anisotropic lattice. In Ref. (2(|, 
we have repeated the improvement procedure of Iwasaki 
on anisotropic lattices, and found that, for small 
anisotropics £ « 1-4, the ^-dependences in the improve- 
ment coefficients are weak, and a sufficient improvement 
is achieved just by fixing the coefficients to Iwasaki's val- 
ues for isotropic lattices, c\ — c\ — c 2 = —0.331. As 
explained in Sec. IIIII we are interested in the case £ = 2. 
Because ^-dependences in the improvement coefficients 
require additional elaborations in numerical simulations, 
such as computation of ^-derivative terms in the equa- 
tion of state (EOS), we fix the improvement coefficients 
to their isotropic values in the following. 



B. Clover quark action on anisotropic lattice 

We employ clover- improved Wilson quarks |27J. On 
anisotropic lattices, the action is given by 



Sf=22 q(x)K(x, y)q(y), 



(4) 



K(x,y) = 6 x>y 

-Kt {(1 - l4 )U 4 (x)S x+4y + (1 + 74)17] (a - i)S x _ ly } 
-k s {(r - li)Ui(x)5 x+ly + (r + 7< )C/?(a: - l)5 x _ ly ) 



-Us \ C t 



^ cr 4l F 4l (x) + rc s ^2 VijFij(x) > 5 x , y . 



(5) 



i>j 



For the field strength F^ v , we use the standard clover- 
leaf definition. Following our previous studies at £ = 1, 
we apply a mean- field improvement for fllf . Ui(x) — » 
Ui{x)/u s and U 4 {x) — > U 4 {x)/ut, where u s and Ut are 
mean links in the spatial and temporal directions. For 
the mean links, we adopt the value estimated from pla- 
quette in one-loop perturbation theory as in our previous 
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FIG. 1: Dispersion relation of free quark on anisotropic lat- 
tices. The energy E(p) — E(0) normalized by £ is plotted as a 
function of spatial momentum at bare quark masses m° q — 0.1, 
0.2 and 0.6 (full, dotted and dashed curves) for anisotropics 
£ = 1, 2 and 4. Larger E(p) — E(0)/£ at Pz/tt ~ 1 correspond 
to smaller £. 



studies. At £ = 2, we obtain 

Wh(ss) = 1 - 1. 154//3, 
Wii(si) = 1 -0.560//?, 



(6) 
(7) 



for spatial and temporal plaquettes. Therefore, we set 
u s = (1 — 1.154//3) 1 / 4 for £ = 2. For the temporal mean- 
field, we adopt ut = 1 because a naively calculated value 
Wuist) 1 / 2 /Wu{ss) 1/4 exceeds 1 for £ >1.6 at our values 
of/?. 

Following Refs. 0, @, we set the spatial Wilson pa- 
rameter to be r = l/£. In this case, the quark dispersion 
relation in physical units preserves the 4-dimensional ro- 
tation symmetry at the tree level, and the tree-level im- 
provement coefficients are free from the terms linear in 
m q UH HE US ■ In a quenched study |3(j , it was shown 
that the fermionic bare anisotropy parameter correspond- 
ing to a fixed £ is well fitted by a quadratic function of 
m q with this choice of r. With dynamical quarks, how- 
ever, terms linear in m q may appear through quark loop 
corrections |29j . 

The clover coefficients, c s and Ct, are functions of r. 
With our choice r = l/£, c s and c* are unity at the tree 
level and their mean-field improved values are given by 



1 



1 



2 ' 

u s uf 



(8) 



We define the bare anisotropy of the fermion field by 

K t U t 



If 



k s u s 



The bare quark mass in units of a s is given by 

1 



/??,, 



7f — 3r. 



(9) 



(10) 



For later convenience, we define n which satisfy the same 
relation with m? q as in the isotropic case: 



- = 2(m 

K V 



4) = 



- 2( 7F + 3r - 4). 



(11) 



We perform chiral extrapolations in terms of 

The relation r = l/£ suggests that spatial doublers 
may appear at large £. The free quark dispersion relation 
for our action is given by 



cosh E(p) = 1 



V + ( 



2~jF 



2(1 



If 



(12) 



where pi = (sinpi)/7^, pi = 2sin(pj/2), and E is in 
units of at while pi is in units of a s 30] . In Fig. ^ we 
plot the energy E(p)/£ for £ = 1, 2 and 4 at m® = 0.1, 
0.2 and 0.6, where £ = 7jt in this approximation. From 
this figure, we expect that doubler effects are weak at 
£ = 2 and m q = (1/k — 1/ k c ) /2 = 0.07-0.8 we study (see 
Sec. ED. 



III. SIMULATION PARAMETERS 

In this paper, we focus on the case of the renormalized 
anisotropy £ = 2. We have shown for finite-temperature 
pure SU(3) gauge theory that this choice of £ is opti- 
mal to reduce scaling violations in equation of state both 
in the high temperature limit and at finite (3; the latter is 
confirmed by a Monte Carlo simulation. It is straightfor- 
ward to analyze the high temperature limit for full QCD. 
We have found that £ = 2 is optimal also with two flavors 
of dynamical quarks and improved glue. 

We perform simulations at (3 = 1.8, 1.9, 2.0, and 2.1 on 
8 3 x 24, 8 3 x 24, 10 3 x 30 and 12 3 x 36 lattices, respectively. 
The lattice spacing is in the range a s ~ 0.28-0.16 fm, and 
hence the spatial lattice size is fixed to be about 2 fm. 
See Sec. lVII for details of the scale determination. At each 
f3, six values of k, corresponding to mps/my ~ 0.6, 0.7, 
0.8, 0.85, 0.9, and 0.92, are simulated. To study lattice 
volume effects, we also perform additional simulations on 
8 3 x 24 and 12 3 x 36 lattices at = 2.0. Our simulation 
parameters are summarized in Table [I] 

We generate gauge configurations by the HMC algo- 
rithm with an even-odd preconditioned BiCGStab quark 
solver |2j. The molecular dynamics time step dt is ad- 
justed to achieve an acceptance rate of about 70-80%. 
Measurements are performed at every 5 trajectories over 
1000-1700 trajectories after 300 thcrmalization trajecto- 
ries, where the length of one trajectory is set to unity. 
Statistical errors of the observables are estimated by the 
jackknife method at each (3 and k with bins of 50 trajec- 
tories. 
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TABLE I: Simulation parameters, "traj." is the number of trajectories used for measurements after 300 thermalization 
trajectories. At f3 = 2.0, besides the main simulation on the 10 3 x 30 lattice, simulations using the same values of k, 7^, 7c 
and the trajectory length are done also on 8 3 x 24 and 12 3 x 36 lattices for a study of finite size effects. 





size 


K 


traj. 


(70, 7fM s ) 












1.8 


8 


x 24 


0.10745 


1700 


(1.70,0.90), 


(1.70,1.10), 


(1 70 1 20) 


(1.75,0.90), 


(1 75 1 10) 


(1 75 1 20) 












(1.85,0.90), 


(1.85,1.10) 
















0.11162 


1700 


(1.70,1.00), 


(1.70,1.20), 


(1.70,1.30), 


(1.75,1.00), 


(1.75,1.20), 


(1.75,1.30) 








0.11582 


1700 


(1.70,1.15), 


(1.70,1.25), 


(1.70,1.30), 


(1.75,1.15), 


(1.75,1.25), 


(1.75,1.30) 








0.12115 


1700 


(1.70,1.25), 


(1.70,1.35), 


(1.70,1.40), 


(1.75,1.25), 


(1.75,1.35), 


(1.75,1.40) 








0.12438 


1700 


(1.70,1.20), 


(1.70,1.30), 


(1.70,1.40), 


(1.70,1.45), 


(1.75,1.30), 


(1.75,1.40), 












(1.75,1.45), 


(1.80,1.20) 
















0.12655 


1700 


(1.70,1.35), 


(1.70,1.40), 


(1.70,1.45), 


(1.75,1.35), 


(1.75,1.40), 


(1.75,1.45) 


1.9 


8 


x 24 


0.1085 


1000 


(1.80,1.00), 


(1.80,1.10), 


(1.80,1.20), 


(1.80,1.30), 


(1.85,1.00), 


(1.85,1.10), 












(1.85,1.20), 


(1.85,1.30), 
















0.1137 


1000 


(1.80,1.15), 


(1.80,1.25), 


(1.80,1.35), 


(1.85,1.15), 


(1.85,1.25), 


(1.85,1.35) 








0.1169 


1000 


(1.75,1.20), 


(1.75,1.30), 


(1.75,1.40), 


(1.80,1.20), 


(1.80,1.30), 


(1.80,1.40), 












(1.85,1.20), 


(1.85,1.30) 
















0.1212 


1000 


(1.75,1.55), 


(1.80,1.25), 


(1.80,1.35), 


(1.80,1.45), 


(1.80,1.55), 


(1.85,1.25), 












(1.85,1.35), 


(1.85,1.45) 
















0.1245 


1500 


(1.70,1.50), 


(1.70,1.60), 


(1.75,1.30), 


(1.80,1.40), 


(1.80,1.50), 


(1.85,1.30), 












(1.85,1.60) 


















0.1260 


1500 


(1.75,1.40), 


(1.75,1.60), 


(1.80,1.50), 


(1.85,1.40), 


(1.85,1.60), 


(1.90,1.50) 


2.0 


10 


x 30 


0.1090 


1000 


(1.80,1.25), 


(1.80,1.35), 


(1.80,1.45), 


(1.85,1.25), 


(1.85,1.35) 






/o3 

(8 


x 24) 


0.1150 


1000 


(1.80,1.45), 


(1.80,1.55), 


(1.85,1.35), 


(1.85,1.45), 


(1.85,1.55), 


(1.95,1.45) 




(12 3 


x 36) 


0.1180 


1000 


(1.80,1.40), 


(1.80,1.50), 


(1.80,1.60), 


(1.85,1.50), 


(1.85,1.60) 










0.1210 


1000 


(1.80,1.45), 


(1.80,1.55), 


(1.80,1.65), 


(1.85,1.45), 


(1.85,1.55), 


(1.95,1.45) 








0.1244 


1500 


(1.70,1.60), 


(1.80,1.50), 


(1.80,1.60), 


(1.80,1.70), 


(1.85,1.55), 


(1.85,1.60), 












(1.90,1.55), 


(1.90,1.60), 


(2.00,1.50) 














0.1252 


1500 


(1.75,1.60), 


(1.75,1.65), 


(1.80,1.60), 


(1.85,1.55), 


(1.85,1.65) 




2.1 


12 3 


x 36 


0.1100 


1000 


(1.80,1.35), 


(1.80,1.55), 


(1.90,1.45), 


(1.95,1.35) 












0.1150 


1000 


(1.80,1.50), 


(1.80,1.60), 


(1.90,1.45), 


(1.90,1.55), 


(1.90,1.65) 










0.1200 


1000 


(1.80,1.65), 


(1.85,1.55), 


(1.90,1.75), 


(1.95,1.50), 


(1.95,1.60) 










0.1225 


1500 


(1.80,1.60), 


(1.80,1.70), 


(1.80,1.80), 


(1.90,1.60), 


(1.90,1.70), (1.90,1.80) 








0.1245 


1500 


(1.80,1.60), 


(1.80,1.80), 


(1.85,1.70), 


(1.90,1.60), 


(1.90,1.70) 





IV. CALIBRATION PROCEDURE 

At each (3 and k, we have to tune the bare anisotropy 
parameters 7_f and 7g such that the renormalized 
anisotropics £p and £g for fcrmionic and gluonic observ- 
ables coincide with each other: 



£f(7f, 7g; P, k) = Cg(7f, 7g; P, k) = C 



(13) 



As discussed in the previous section, we study the case 
£ = 2. For this purpose, we measure £p and £g at sev- 
eral values of (7_f,7g) at fixed k and /3, and determine 
the point where Eq. (|13(l is satisfied by an interpolation 
in 7i? and 70- Let us denote the resulting values of 
and 7g for £ = 2 as j f ((3,k) and 7 G (/3, k). Finally, we 
parametrize 7^ and 7 G as functions of (3 and k for use 
in future studies of heavy quark systems and thermody- 
namics of QCD. 



We measure £g by Klassen's method [Tfj|: 

R s (x,y) = R t (x,&y), (14) 

where 

W ss (x,y) 



R s {x,y) 
Rt(x,t) 



W ss (x + l,y) 

Wst(x,t) 
W st (x + l,t)' 



(15) 
(16) 



are the ratios of spatial-spatial and spatial-temporal Wil- 
son loops, W ss (x,y) and W s t(x, t), respectively. We de- 
termine £g by minimizing 



x.y 



(R s (x,y) ~ Rt{x,£ G y))< 

(AR s y + (AR t y 



(17) 
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with AR S and AR t the statistical errors of R s and Rt- 
To avoid short range lattice artifacts, x and y should not 
be too small. Practical range of x and y will be discussed 
later. 

For £p we use the relativistic dispersion relation of 
mesons: 



,2 i P 



E(p) 2 = m 2 + ^+0(f), 

SF 



(18) 



where E and m are the energy and mass in units of at, 
and p = 2irn/L s , with L s the spatial lattice size, is the 
spatial momentum in units of a s . We evaluate E and m 
from a cosh fit of the meson two-point correlation func- 
tion, 



C(p,t) = £<O(x,i)Ot(0,0)e^> 



(19) 



0(x,t) = ^2cf>(f)(l)'(z)q(x + y,t)rq(x + z,t) (20) 
yz 

In this paper, we study pseudoscalar (PS) and vector (V) 
mesons consisting of sea quarks only: T = 75 for PS and 
r = 7j for V. Quark fields are smeared by a function 
4>{x) to enhance ground state signals at short distances. 
For the "smeared" quark field, we adopt an exponential 
smearing function of form 

<f>(x) = aexp(-p|x|) fori ^ 0, 0(0) = 1, (21) 

where the coefficients a and p are adopted from a previous 
study 0]. The "point" quark field corresponds to <p(x) = 
5g g. In our calculation of meson two-point function, the 
sink operator is always the point-point type, while, for 
the source operator, we study point-point, point-smeared 
and smeared-smeared cases. We find that the smear- 
smear source operator leads to the earliest plateau with 
small errors. Therefore, we adopt the smear-smear source 
operator. 

In principle, we may adopt different observables to de- 
fine the renormalizcd anisotropics. Away from the con- 
tinuum limit, different choices will lead to 0(a) differ- 
ences in the calibration results. To study this problem, 
we compare the calibration results using £p from PS and 
V meson dispersion relations. We denote these results 
for the calibrated bare anisotropies as ( , y F (PS),j G (PS)) 
and {j F (V), J G (V)), respectively. In Sec. IV CI we show 
that they tend to converge together toward the contin- 
uum limit. In future applications of the present work, 
different sets of (7^,7^) wm be useful for estimating 
systematic errors due to the continuum extrapolation, in 
complicated physical observables, such as the equation of 
state. 

In a previous study of quenched QCD the ratio 
of temporal and screening masses of the PS meson was 
used to determine £p. We study the difference between 
our procedure and the mass ratio method in Appendix 
A. We find that both methods give consistent values of 
£f when the quarks are not too heavy [mps/mv ^0.75 
(0.8) at P >2.0 (2.1)]. 
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FIG. 2: A typical for a determination of ^ at /3 = 2.0, 
k = 0.1244 and (-jg,1fu s ) = (1.85, 1.60). The £g shows the 
minimizing point of L(£q) defined by Eq |17l Different sym- 
bols represent the results obtained on 8 3 x 24, 10 3 x 30 and 
12 2 x 36 lattices. For filled/open symbols, £(£g) is evaluated 
with/without the y = 1 data. 



V. CALIBRATION RESULTS 
A. £g from matching of Wilson loop ratios 

We determine the renormalized gauge anisotropy £q by 
minimizing the function £(£g) defined by Eq. i|17|) . We 
interpolate Rt(x,t) by a cubic spline in terms of t. To 
remove short range lattice artifacts, we evaluate £(£g) 
with x and y which satisfy x x y > M and examine 
the M dependence. The upper limit on x and y is set 
by requiring that the statistical error does not exceed 
the central value for the Wilson loop ratio. Varying the 
upper limit hardly changes results for £g ■ Filled symbols 
in Fig. 13 show typical results of £g as a function of M = 
min(a; x y). We find that, at this simulation point, £g is 
reasonably stable when x x y is larger than about 4 . 

Since the condition xxy > M does not exclude small x 
or y, which can be an additional origin of short distance 
effects, we study whether £g are affected by small values 
of x or y by removing them. The results of £g using L(£g) 
without data at y = 1 are plotted with open symbols in 
Fig. |21 We find that, although clear deviations from the 
filled symbols are observed at small M, the effects of 
y = 1 data are within 1% at M>4 where £g becomes 
stable. Therefore, placing a condition on min(x x y) is 
sufficient to obtain a stable value for £q. Similar results 
are obtained at other simulation points. 

Results obtained on different lattice volumes (8 3 x 24, 
10 3 x 30 and 12 2 x 36) are also shown in Fig.|2| With our 
lattices, no finite volume effects are visible in the values 

From these studies, we adopt min(a; x y) = 3, 3, 4, 
and 5 at f3 = 1.8, 1.9, 2.0 and 2.1, respectively, in the 
subsequent analyses. 
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FIG. 3: Effective mass of meson states with various momenta 
obtained at p = 2.0, k = 0.1244, (■y G ,-yFU s ) = (1.85, 1.60) on 
the 10 3 x 30 lattice. The left and right panels are the results 
for pseudoscalar and vector mesons. 
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FIG. 5: Volume dependence of a mesonic dispersion relations 
at f3 = 2.0, k = 0.1244 and (-ya^FU,) = (1.85, 1.60). Filled 
and open symbols show the results on 12 3 x 36 and 8 3 x 24 
lattices with PS and V channels. Other condition is same as 
Fig0]and fit results on 10 3 x 30, which are the same as FiglH 
are shown with dotted lines. 



Using data at ft — (0,0,0), (1,0,0) and its permu- 
tations, we fit E(p) with the leading formula E(p) 2 = 
m 2 +p^/£,p to determine £f- The fits are shown by dot- 
ted lines in Fig. 0] In Fig. [SJ data obtained on 8 3 x 24 
and 12 2 x 36 lattices are compared with the fit results 
on the 10 3 x 30 lattice. We find that the data are well 
explained by the fit results. Indeed, the slope obtained 
for the three lattice sizes are consistent: £f = 2.044(59), 
2.020(44) and 2.012(39) on 8 3 x 24, 10 3 x 30 and 12 2 x 36 
lattices, respectively. This confirms that the spatial lat- 
tice size >1.6 fm is sufficiently large to suppress finite 
volume effects in £f in the range of quark masses we 
study. 



FIG. 4: Sample results for the dispersion relation of pseu- 
doscalar and vector mesons at /3 — 2.0, k = 0.1244 and 
(iGj'yFUs) = (1.85,1.60). Dotted lines show fit results from 
ft = (0,0,0) and (1,0,0). 



B. £f from meson dispersion relations 

We determine the renormalized quark anisotropy £f 
from the meson dispersion relation. We calculate the 
meson energy E(p) at the spatial momenta p = 2irn/L s 
with ft = (0, 0,0), (1,0,0), (1,1,0), (1,1,1), and their per- 
mutations. In Fig. |2| we plot typical data for effective 
energy defined by 

C(p,t) = cosh[E eS (p,t)(N t /2-t)] 
C&t + l) cosh[E eS (p,t){N t /2-t-l)Y 1 ' 

obtained from the smear-smear correlators. Typical re- 
sults for the energy E(p) are shown in Fig. 0] 



C. Bare anisotropics at £ = 2 (7J and 7^) 

Figures show typical results for bare anisotropies 
at each (3, obtained at the third and fifth heaviest quark 
masses (mps/mv ~ 0.85 and 0.70). We find that, for 
the range of parameters we study, we can fit the data 
assuming an ansatz linear in jp and jg : 

£f = a F + bpj F + c F 7 G , (23) 
= a G + b G j F + cgIg- (24) 

Results of the least \ 2 fits are also shown in Figs. 

From the condition £f(7f,7 g ) = £g(7f>7 g ) = 2 > we 
obtain 7^ and 7 G for £ = 2 as functions of the coefficients, 
ajy, • • • , cq. We determine their errors using the error 
propagation formula where the errors for ap, • ■ • , cq are 
estimated from the error matrix of the least \ 2 fits for 
£f and £g- The results are summarized in Table ITT1 

To confirm the magnitude of errors, we study 
X 2 totai/ N df = (x F :2 + X G - : 2)/ 2N df, as a function of 7f and 
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TABLE II: Bare anisotropy parameters calibrated to £ = 2. 







7&(PS) 


2 / 7Vr 

X /N df 


7f(PS) 


2 / A T 

X /N df 


ib (v) 


2 / a r 

X /N df 


7f(V) 


2 / AT 

X /Ndf 


1.8 


0.10745 


1.742(11) 


12.4/ 5 


1.499(12) 


6.7/ 5 


1.7748(93) 


12.4/ 5 


1.2783(82) 


3.1/ 5 




0.11162 


1.734(12) 


2.7/ 3 


1.571(17) 


2.6/ 3 


1.740(18) 


2.7/ 3 


1.298(21) 


2.0/ 3 




0.11582 


1.725(11) 


2.2/ 3 


1.682(20) 


1.2/ 3 


1.768(20) 


2.2/ 3 


1.464(43) 


2.4/ 3 




0.12115 


1.7282(81) 


0.9/ 3 


1.784(11) 


1.4/ 3 


1.757(15) 


0.9/ 3 


1.613(14) 


7.8/ 3 




0.12438 


1.703(13) 


2.1/ 5 


1.879(16) 


3.3/ 5 


1.769(14) 


2.1/ 5 


1.676(28) 


2.5/ 5 




0.12655 


1.708(22) 


1.1/ 3 


1.895(15) 


4.9/ 3 


1.760(22) 


1.1/ 3 


1.786(13) 


1.4/ 3 


1.9 


0.10850 


1.8068(80) 


1.6/ 5 


1.554(12) 


5.1/ 5 


1.8093(56) 


1.6/ 5 


1.387(15) 


5.8/ 5 




0.11370 


1.7971(94) 


2.9/ 3 


1.675(12) 


1.3/ 3 


1.8126(60) 


2.9/ 3 


1.5615(84) 


1.6/ 3 




0.11690 


1.760(11) 


1.3/ 5 


l. 753(13) 


7.2/ 5 


1.7945(74) 


1.3/ 5 


1.6085(82) 


3.8/ 5 




0.12120 


1.773(10) 


4.0/ 5 


1.858(13) 


1.6/ 5 


1.8013(75) 


4.0/ 5 


1.742(13) 


7.3/ 5 




0.12450 


1.7631(72) 


2.1/ 4 


1.927(13) 


2.2/ 4 


1.7908(85) 


2.1/ 4 


1.807(30) 


1.5/ 4 




0.12600 


1.7629(93) 


2.8/ 3 


1.987(17) 


1.8/ 3 


1.7811(77) 


2.8/ 3 


1.918(12) 


4.9/ 3 


2.0 


0.10900 


1.8243(40) 


2.9/ 2 


1.6931(70) 


4.1/ 2 


1.8319(38) 


2.9/ 2 


1.6266(95) 


1.3/ 2 




0.11500 


1.8288(52) 


1.5/ 3 


1.8015(85) 


0.9/ 3 


1.8331(76) 


1.5/ 3 


1.7190(98) 


6.7/ 3 




0.11800 


1.8243(65) 


5.9/ 2 


1.8907(74) 


0.2/ 2 


1.8315(88) 


5.9/ 2 


1.8268(83) 


0.1/ 2 




0.12100 


1.8225(73) 


7.5/ 3 


1.935(19) 


3.3/ 3 


1.8278(71) 


7.5/ 3 


1.903(25) 


2.5/ 3 




0.12440 


1.8120(73) 


12.1/ 6 


2.002(19) 


1.7/ 6 


1.8299(92) 


12.1/ 6 


1.938(27) 


3.1/ 6 




0.12520 


1.8169(76) 


2.2/ 2 


2.026(14) 


0.4/ 2 


1.8289(78) 


2.2/ 2 


1.949(58) 


5.5/ 2 


2.1 


0.11000 


1.8814(88) 


2.7/ 1 


1.796(10) 


1.4/ 1 


1.8827(79) 


2.7/ 1 


1.760(11) 


0.4/ 1 




0.11500 


1.8678(70) 


4.5/ 2 


1.8932(74) 


17.5/ 2 


1.8722(73) 


4.5/ 2 


1.8501(85) 


19.4/ 2 




0.12000 


1.8673(87) 


0.0/ 2 


1.970(18) 


1.2/ 2 


1.871(16) 


0.0/ 2 


1.954(49) 


0.6/ 2 




0.12250 


1.8559(58) 


2.6/ 3 


2.032(23) 


3.6/ 3 


1.8603(62) 


2.6/ 3 


2.004(20) 


6.0/ 3 




0.12450 


1.8517(55) 


1.6/ 2 


2.043(14) 


1.0/ 2 


1.8615(68) 


1.6/ 2 


1.980(29) 


4.8/ 2 




B=1.8,k=0.12438 
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.70 




■ Yc- 


.75 




Ay,- 


.80 ~ 


' ' ' i ' ' ' ' i ' 

- p=1.8,K=0.12438 


I I 




X Z /N dl =2.1/5 












A - 


•v 






.70 




■7c- 


.75 




Av 


.80 



1.3 



2.2 

2.1 

2.0 

| 1.9 

1.8 

1.7 

1.6 
2.2 

2.1 

„ 2.0 

1.9 

1.8 

1.7 



P=1.9,k=0.1169 




I X Z /N al =7.2/5 






• T.-175 




■ 7g=.1.80 






I I I I ' ' ' ' I 

- P=1.9,k=0.1169 


I I I I I I I I 


X Z /N„,=1.3/5 










•V.-1-75 




■ y G =l.ao _ 




AV.-1-85 



2. 
2. 
2 

ro 2. 

Q_ 

i < 

ULP 1 ■' 
1.1 

1.: 



1.3 



2.1 
, 20 
1.9 
1.8 
1.7 



P=1.9,k=0.1245 




Z Z /N al =2.2/4 




I i i i i I i i i i I 


• 7.-1.70 - 
■ 7.-1-75 

AYg-1-80 : 

OTg-1-85 


P=1.9,k=0.1245 


H l | 


X 2 /N d ,=2.1/4 




. « ' A 

• . A .. • 

. - * " " • 




•Y„-1.70 
■ V.-1.75 : 

AV G =1.80 

07.-1.85 



1.4 1.5 



FIG. 6: £g and £f as functions of (7g,7fUs) at /3 = 1.8 
for the third and fifth heaviest k. For £f, results from the 
pseudoscalar dispersion relation are shown. Lines represent 
the results of the fits and Elt . 



7 G , where Xp/ G;2 is the \ 2 value for a fit of £ F / G ("f F , j' G ) 
data to £f/g = 2 + 6(7^? — 7f) + C (7 G ~ 7g) f° r given val- 
ues of (7F)7g)- This quantity measures to which extent 
£ = 2 is achieved by Wilson loops and the meson corre- 
lation function at (7f,7g)- The minimum of Xtotai/^df 



FIG. 7: The same as Fig. ©but at /3 = 1.9. 



locates at (7^ , 7 G ). A typical result is plotted in Fig. 1101 
We find that errors estimated from a unit increase of 
Xtotai/Ndf are consistent with those listed in Table ITT1 

In later applications, it will be convenient to 
parametrize 7^ and 7 G as functions of (3 and K. Fig- 
ures ^2 and E| show the parameter dependence of 7^ 
and 7 G . We adopt the general quadratic ansatz in (3 and 
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FIG. 8: The same as Fig. ©but at f3 = 2.0. 
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FIG. 10: xiotai/Ndf for £ = 2 fits at /3 = 2.0 and k = 
0.121. Results using £f(PS) for the fermionic anisotropy are 
shown. The minimum of Xtotai/Ndf is 1.348 at (u s 7f,7g) = 
(m s 7f)7g) = (1.5605, 1.8225). Curves on the base plane is a 
contour map of Xtotai/Ndf- The minimum point is marked by 
"x" on the contour map. [Note that Ndf is larger than that 
for the fits 12311 and 12411 summarized in Table ITU because a F 
and ac are not free in this calculation.] 



FIG. 9: The same as Fig. ©but at = 2.1. 




K, 



1.6 



7 * = A F + B F p' + C F (3' 2 

+D f /3'k' + EpK 1 + F f k' 2 , 
7c, = A G + B G f3> + C G {3' 2 

+D g 0'k' + E g k' + F G n' 2 , (26) 



(25) 



where 0' = /3 - 2.0 and k' = K - 0.12. For 7^ (PS) and 
7^ (PS), we find 

A F = 1.9097(41), B F = 0.746(38), 

C F = 0.04(25), Dp = -13.9(4.1), 

E F = 20.14(75), Fp = -1.(83.), (27) 

A G = 1.8210(28), B G = 0.435(22), 

C G = -0.24(17), D G = 3.2(2.9), 

E G = -1.69(44), F G = -69.(57.), (28) 



with X 2 /N df = 38.0/17 and X 2 /N df 
7 >(V) and 7 £(V), 



19.2/17, and, for 
1.8434(48), Bp = 1.204(46), 
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FIG. 11: 7g(PS) and 7g(V) corresponding to £ = 2. Curves 
are the results of the global fit I26H with parameters 1281 and 
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FIG. 12: 7]? (PS) and 7f(V) corresponding to £ = 2. Curves 
are the results of the global fit I25H with parameters 127H and 
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FIG. 13: 7^(PS) vs quark mass at various /3. Lines are guide 
of the eyes. 
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FIG. 14: Relative differences between 7^(PS) and 7f(V), and 
between 7g(PS) and 7g(V), as functions of /3 and 1/k. 



C f = -0.93(29), D F = -21.9(4.7), 

E F = 25.94(97), F F = 488.(92.), (29) 

Aq = 1.8311(29), Bq = 0.348(22), 

C G = 0.14(20), Dq = 0.5(2.9), 

E G = -1.19(47), Fq = -61.(58.), (30) 

with x 2 /N df = 76.0/17 and x 2 /N df = 14.6/17, using the 
values of J F i G and their errors listed in Table. |HJ The 

errors for the coefficients are estimated from the x 2 error 
matrix. These fits are shown in Figs.llllandll2lbv dotted 
lines. 

In Fig. 1131 we plot 7^ as a function of the dimcnsionlcss 
quark mass m q = (1/k — 1/k c )/2 using k c determined in 
Sec. EU Although the range of m q is not quite close to 
the chiral limit, our values of 7^ suggest a strong linear 
dependence in m q . This result is in clear contrast to the 
case of quenched QCD in which 7^ is well fitted by a 
quadratic ansatz in m q motivated from a tree level ex- 
pression of 7^ |3C|. As mentioned in Sec. Ill Bl we expect 
linear corrections from higher order quark loops even if 
the linear terms are removed at tree level [25^ . Our result 
provides us with an example which confirms this expec- 
tation. 

Finally, we study the 0(a) differences among the cal- 
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ibration results using PS and V mesons. We plot rela- 
tive differences between 7^ (PS) and 7p(V), and between 
7g(PS) and J G (V), in Fig. El as functions of (3 and 1/k. 
Errors are estimated neglecting the correlation between 
PS and V determinations. We find that the differences 
tend to vanish as (3 is increased. At (3 > 2.0, the differ- 
ences are less than 5% for 7^ and 1% for -f G . 

VI. PHYSICAL QUANTITIES AT £ = 2 

In this section, we interpolate the measurement results 
to the calibration points corresponding to £ = 2 to esti- 
mate the scale and several other basic properties of our 
lattice. We also test the effects of the two calibration 
results using PS and V meson dispersion relations on the 
continuum extrapolation of physical quantities. 

For the interpolations to (7^,7^) at each (/?, k), we 
adopt a linear ansatz 

y — a + bjo + crfF (31) 

when the range of 7^ is less than 0.3. When max(7i?) — 
mm(7F ) > 0.3, we adopt a quadratic ansatz 

y = a + bj G + crf F + d-f 2 Fl (32) 

because, in this case, the linear ansatz sometimes fails 
to explain the data (x 2 /N df ~ O(10)-O(100)). We find 
that terms quadratic in 7^ do not improve the fits. We 
confirm that this quadratic ansatz leads to a result con- 
sistent with the linear ansatz if max(7i?) — min(7i?) < 0.3. 

Several physical quantities thus interpolated to £ = 
2 are summarized in Table. IIIII where the results from 
the quadratic ansatz (|32|) are marked by "*" on k. The 
errors are estimated by quadratically averaging over the 
contributions from the x 2 error matrix for the fits (|31[) or 
(|32fl , and from the errors for 7^ and j G . Results adopting 
two alternative choices for the £ = 2 point — 7^y G (PS) 
from the pseudo-scalar dispersion relation and 7jty G (V) 
from the vector dispersion relation — are labeled by (PS) 
and (V). 
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FIG. 15: Plaquettes at £ = 2 and 1. Different points at the 
same /3 with the same symbol represent the results obtained at 
different k (larger Wu(ss) and Wn(st) correspond to larger 
k). Dashed lines are the results of 1-loop perturbation theory. 
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FIG. 16: m 2 PS at 7j5-/g(P>3) for £ = 2 as a function of 1/k. 
Full and dotted curves show quadratic and linear fits using 
lightest four and three data points, respectively. 



A. Plaquette 

Figure El shows the plaquette expectation values 
W n (ss) and W n {st) at £ = 2 adopting 7* /G (PS). Re- 
sults adopting 7£yg(V) are similar. As a reference point, 
plaquette values on an isotropic lattice using the same 
action are also plotted 0]. Different points at the same 
[3 with the same symbol represent the results obtained at 
different k. 

Numerical results of plaquettes are compared with 
their 1-loop values at £ = 2, Eqs. © and (J7J, shown 
by dashed lines in the Figure. For (3 we adopt the bare 
value. We find that, as in the case of isotropic lattices, 
the plaquettes agree with the perturbative calculation 
within about 10% at this value of (3. This confirms our 
choice of the clover coefficients Ct and c s to this accuracy. 



B. Light meson spectrum and the lattice scale 

In Fig. 1161 we summarize the values of rn? PS at £ = 2 
adopting 7^y G (PS), as a function of 1/k. Carrying out 
chiral extrapolations in which the lightest four points are 
fitted to a quadratic ansatz, we obtain the chiral point 
k c (PS) listed in Table llVl The first errors are statistical, 
while the second ones are systematic errors estimated 
from the difference with the results of linear fits to the 
lightest three points. The values of k c (V) are obtained 
similarly, adopting 7^y G (V) as the point for £ = 2. We 
find that k c (PS) and k c (V) are consistent with each other 
within the present statistical accuracy. 

We determine the scale of our lattices from the p 
meson mass m p = 771.1 MeV at the physical point 
mps/mv — fn^/mp — 135.0/771.1. FigurelT7lshows mi/ 
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TABLE III: Physical quantities interpolated to £ = 2. Results marked by "*" are from the quadratic ansatz 13211 . while other 
results are from the linear ansatz <|31^ . Results of two alternatives for the £ = 2 point, 'yp/aiPS) and 7wg(V), are labeled by 
(PS) and (V), respectively. 





K 


m PS (PS) 


X 2 /N df 


m y (PS) 


X 2 /N df 


m PS /m v (PS) 


X 2 /N df 


r /a s (PS) 


X 2 /N df 






mp S (V) 


X 2 /N df 


m v (V) 


X 2 /N df 


■m PS /mv(V) 


X 2 /N df 


r /a s (V) 


X 2 /N df 


1.8 


0.10745* 


1.4395(66) 


8.2/ 4 


1.5760(75) 


2.6/ 4 


0.91339(37) 


1.3/ 4 


1.299(15) 


2.6/ 4 






1.5601(55) 


8.2/ 4 


1.7178(64) 


2.6/ 4 


0.90828(43) 


1.3/ 4 


1.320(19) 


2.6/ 4 




0.11162* 


1.2431(62) 


1.1/ 2 


1.3985(75) 


3.1/ 2 


0.88900(57) 


2.8/ 2 


1.346(18) 


1.5/ 2 






1.3009(39) 


1.1/ 2 


1.4720(45) 


3.1/ 2 


0.88382(65) 


2.8/ 2 


1.355(19) 


1.5/ 2 




0.11582 


1.0302(57) 


7.9/ 3 


1.2043(78) 


7.1/ 3 


0.8556(10) 


3.1/ 3 


1.413(11) 


3.0/ 3 






1.079(12) 


7.9/ 3 


1.275(17) 


7.1/ 3 


0.8458(23) 


3.1/ 3 


1.406(19) 


3.0/ 3 




0.12115 


0.7635(15) 


2.9/ 3 


0.9679(25) 


5.5/ 3 


0.7888(12) 


3.0/ 3 


1.5210(64) 


8.5/ 3 






0.7660(26) 


2.9/ 3 


0.9906(42) 


5.5/ 3 


0.7732(19) 


3.0/ 3 


1.557(12) 


8.5/ 3 




0.12438* 


0.5880(25) 


5.7/ 4 


0.8165(36) 


7.4/ 4 


0.7200(22) 


6.1/ 4 


1.605(10) 


8.1/ 4 






0.5448(73) 


5.7/ 4 


0.7973(46) 


7.4/ 4 


0.6834(74) 


6.1/ 4 


1.701(13) 


8.1/ 4 




0.12655 


0.4347(71) 


9.1/ 3 


0.6841(45) 


4.3/ 3 


0.6356(80) 


3.9/ 3 


1.730(11) 


4.2/ 3 






0.3838(62) 


9.1/ 3 


0.6627(46) 


4.3/ 3 


0.5805(70) 


3.9/ 3 


1.796(10) 


4.2/ 3 


1.9 


0.10850* 


1.3255(56) 


4.5/ 4 


1.4428(65) 


2.3/ 4 


0.91851(50) 


1.5/ 4 


1.5015(66) 


18.2/ 4 






1.4045(79) 


4.5/ 4 


1.5330(90) 


2.3/ 4 


0.91599(54) 


1.5/ 4 


1.5101(73) 


18.2/ 4 




0.11370 


1.0672(38) 


7.5/ 3 


1.2002(46) 


5.7/ 3 


0.88860(95) 


10.3/ 3 


1.6041(80) 


2.4/ 3 






1.0977(27) 


7.5/ 3 


1.2393(32) 


5.7/ 3 


0.88535(59) 


10.3/ 3 


1.6140(49) 


2.4/ 3 




0.11690 


0.9085(30) 


4.9/ 5 


1.0557(42) 


8.9/ 5 


0.86105(98) 


9.4/ 5 


1.6376(71) 


4.3/ 5 






0.9304(20) 


4.9/ 5 


1.0891(27) 


8.9/ 5 


0.85428(61) 


9.4/ 5 


1.6742(44) 


4.3/ 5 




0.12120* 


0.6825(15) 


5.8/ 4 


0.8460(27) 


5.7/ 4 


0.8066(17) 


2.1/ 4 


1.795(10) 


1.8/ 4 






0.6793(13) 


5.8/ 4 


0.8516(25) 


5.7/ 4 


0.7976(18) 


2.1/ 4 


1.821(10) 


1.8/ 4 




0.12450* 


0.4859(24) 


5.9/ 3 


0.6720(21) 


0.6/ 3 


0.7234(28) 


2.7/ 3 


1.9708(70) 


19.1/ 3 






0.4593(73) 


5.9/ 3 


0.6599(40) 


0.6/ 3 


0.6958(81) 


2.7/ 3 


2.016(18) 


19.1/ 3 




0.12600 


0.3804(61) 


33.6/ 3 


0.5839(39) 


8.5/ 3 


0.6517(79) 


10.7/ 3 


2.119(14) 


1.6/ 3 






0.3559(43) 


33.6/ 3 


0.5722(30) 


8.5/ 3 


0.6205(56) 


10.7/ 3 


2.165(10) 


1.6/ 3 


2.0 


0.10900 


1.1913(26) 


0.9/ 2 


1.2830(30) 


1.7/ 2 


0.92857(36) 


0.6/ 2 


1.7792(53) 


3.9/ 2 






1.2143(35) 


0.9/ 2 


1.3097(40) 


1.7/ 2 


0.92725(39) 


0.6/ 2 


1.7844(55) 


3.9/ 2 




0.11500 


0.9057(18) 


2.2/ 3 


1.0131(24) 


6.2/ 3 


0.89420(51) 


8.3/ 3 


1.9326(45) 


8.2/ 3 






0.9201(22) 


2.2/ 3 


1.0323(30) 


6.2/ 3 


0.89158(69) 


8.3/ 3 


1.9471(66) 


8.2/ 3 




0.11800 


0.75539(97) 


1.2/ 2 


0.8696(14) 


0.3/ 2 


0.86874(67) 


0.4/ 2 


2.0380(48) 


5.1/ 2 






0.7604(13) 


1.2/ 2 


0.8778(18) 


0.3/ 2 


0.86623(92) 


0.4/ 2 


2.0553(68) 


5.1/ 2 




0.12100 


0.59609(94) 


26.3/ 3 


0.7225(14) 


14.6/ 3 


0.8251(15) 


0.6/ 3 


2.229(11) 


1.5/ 3 






0.59505(95) 


26.3/ 3 


0.7231(15) 


14.6/ 3 


0.8229(19) 


0.6/ 3 


2.241(11) 


1.5/ 3 




0.12440 


0.3912(35) 


11.5/ 6 


0.5397(23) 


20.9/ 6 


0.7254(41) 


11.6/ 6 


2.527(13) 


14.6/ 6 






0.3784(50) 


11.5/ 6 


0.5319(31) 


20.9/ 6 


0.7116(57) 


11.6/ 6 


2.570(18) 


14.6/ 6 




0.12520 


0.3355(37) 


8.5/ 2 


0.4906(28) 


5.6/ 2 


0.6843(56) 


0.1/ 2 


2.599(16) 


0.3/ 2 






0.315(15) 


8.5/ 2 


0.4819(67) 


5.6/ 2 


0.655(23) 


0.1/ 2 


2.660(47) 


0.3/ 2 


2.1 


0.11000 


1.0442(32) 


3.2/ 1 


1.1150(36) 


0.2/ 1 


0.93648(52) 


3.4/ 1 


2.2324(72) 


13.1/ 1 






1.0545(32) 


3.2/ 1 


1.1264(36) 


0.2/ 1 


0.93610(47) 


3.4/ 1 


2.2372(61) 


13.1/ 1 




0.11500 


0.8108(13) 


9.0/ 2 


0.8898(16) 


4.8/ 2 


0.91143(49) 


0.0/ 2 


2.4000(63) 


1.5/ 2 






0.8162(14) 


9.0/ 2 


0.8963(17) 


4.8/ 2 


0.91097(59) 


0.0/ 2 


2.4134(68) 


1.5/ 2 




0.12000 


0.56448(90) 


1.3/ 2 


0.6577(14) 


1.3/ 2 


0.8584(11) 


0.5/ 2 


2.6992(94) 


0.3/ 2 






0.5644(14) 


1.3/ 2 


0.6582(24) 


1.3/ 2 


0.8577(20) 


0.5/ 2 


2.705(18) 


0.3/ 2 




0.12250 


0.4313(14) 


5.8/ 3 


0.5339(12) 


3.6/ 3 


0.8075(28) 


5.3/ 3 


2.875(10) 


10.0/ 3 






0.4294(13) 


5.8/ 3 


0.5335(13) 


3.6/ 3 


0.8044(25) 


5.3/ 3 


2.881(11) 


10.0/ 3 




0.12450 


0.3015(23) 


10.4/ 2 


0.4220(19) 


11.9/ 2 


0.7134(32) 


7.9/ 2 


3.177(15) 


1.6/ 2 






0.2906(48) 


10.4/ 2 


0.4149(34) 


11.9/ 2 


0.6997(63) 


7.9/ 2 


3.229(27) 


1.6/ 2 
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FIG. 17: m v as a function of m 2 PS for £ = 2 at 7f/g(PS) 
(filled symbols) and at 7wq(V) (open symbols). The left- 
most symbols are the results of quadratic extrapolations to 
the physical point. The lines are guide for the eyes. 



TABLE IV: Critical hopping parameter k c . The numbers in 
the first bracket are statistical errors, and those in the second 
bracket are systematic errors from the chiral extrapolation. 



p 


l/«c(PS) 


1/MV) 


1.8 


7.708(43)(±o 3 ) 


7.763(80) (+o 5 ) 


1.9 


7.775(39) (If) 


7.789(46) (±|°) 


2.0 


7.836(26)(±o 5 ) 


7.862(50)(± 2 , 2 ) 


2.1 


7.888(14)(tf 7 ) 


7.896(24) (to 11 ) 



as a function of m 2 PS at 7^ G (PS) and 7^/ G (V). We find 
that the masses at 7^ G (PS) and 7^/ G (V) are slightly 
different on coarse lattices. The difference rapidly de- 
creases with increasing (3. 

Extrapolation to the physical point is done adopting a 
quadratic ansatz to the lightest four data points at each 
P. We note that, with the present statistics, two results 
of mv at 7^y G (PS) and 7^y G (V), extrapolated to the 
physical point, are consistent with each other already on 
the coarsest lattice. High statistics simulations directly 
at 7^ G (PS) and 7^y G (V) may resolve the difference at 
the physical point. 

The resulting lattice scale is summarized in Table Ivl 



TABLE V: Lattice scale determined from m p at the physical 



point. 


13 


a s (PS) [GeV" 1 ] 


a fl (V) [GeV" 1 ] 


L s a s [fm] 


1.8 
1.9 
2.0 
2.1 


1.395(28)(^ 5 ) 
1.185(21)(t 8 °) 

0.957(15) (to 18 ) 
0.824(10)(lg 2 ) 


1.408(26) (iS 7 ) 
1.178(17)(±^ 5 ) 
0.986(26) (l® 1 ) 
0.838(17)(15°) 


2.2 
1.9 
1.9 
2.0 



TABLE VI: Parameters for the calculation of the static quark 
potential. 



p 




K 


Wsmcar 


fit-rang 


;e in t 


fit-ran; 


je in 


1.8 


0.10745 


- 0.11582 


l 


3 - 


6 


V2- 


2V3 




0.12115 


- 0.12655 


l 


3 - 


6 


V2- 


3V2 


1.9 


0.1085 


- 0.1169 


2 


3 - 


6 


V2- 


3^2 




0.1212 


- 0.1260 


2 


4 - 


7 


V2- 


3V3 


2.0 


0.1090 


- 0.1180 


3 


4 - 


7 


V2~ 


3V3 




0.1210 


- 0.1252 


3 


5 - 


9 


V3- 


- 6.0 


2.1 


0.1100 


- 0.1200 


4 


5 - 


9 


V3- 


- 6.0 




0.1225 


- 0.1245 


4 


6 - 


11 


2.0 - 


3^6 



3.5 




1.0 1 — 

0.0 0.5 1.0 1.5 2.0 

m ps 2 

FIG. 18: Sommer scale at £ = 2 as a function of the PS 
meson mass. The leftmost symbols are the results of quadratic 
extrapolations to the physical point. The lines are guide for 
the eyes. 

Systematic errors (second errors) in the table are esti- 
mated from a comparison with linear fits using the light- 
est three points. 



C. Static quark potential and Sommer scale 

We extract the static quark potential V(r) from a fit 
W{r, t) = C{r) exp [-V{r)t] , (33) 



TABLE VII: Sommer scale at £ = 2, extrapolated to the 
physical point. 



J3 rp/aJPS) r /a s (V) 

1.8 1.843(21)f+° 7 ) 1.892(17)(+° 6 ) 

1.9 2.269(23) (±° rr ) 2.331(18)(tg 6 ) 

2.0 2.818(27) (±£ 2 ) 2.870(41) (±? 8 ) 

2.1 3.367(23) (±^) 3.377(39) (±? 18 ) 
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1.5 



is given by 



1.3 



E 1-2 - 



1.0 - 



•m 


s /m„=0.175 


PS 


H 


0.60 


PS 


A 


0.70 


PS 


Om 


s /m v =0.175 


V 


□ 


0.60 


V 


A 


0.70 


V 



0.9 



0.2 



0.4 



0.6 

m„ 



1.0 



FIG. 19: Sommer scale at mps/my = 0.7, 0.6 and 0.175 (the 
physical point) as a function of the lattice spacing. Errors are 
statistical. 



1.10 
1.05 
1.00 
0.95 
0.90 



physical point 
•WPS) 

ot™00 



ni_ 

ft.s 



'1.65 + a£ 



(36) 



o.o 



0.2 



0.4 



0.6 



Results of ro/a s interpolated to £ = 2 are listed in Ta- 
bic [ffi] and shown in Fig.^3] Extrapolating to the physi- 
cal point, we obtain the values summarized in Table IVTT1 
where the central values are from quadratic fits in mp s 
using the lightest four k's and the systematic errors are 
estimated from the difference with a linear fit using the 
lightest three k's. 

The Sommer scale at mp S /m v — 0.7, 0.6 and 0.175 
(the physical point) from a quadratic fit is plotted in 
Fig-EH as a function of the lattice spacing. We find that 
the difference between the calibrations using the PS and 
V meson dispersion relations becomes smaller toward the 
continuum limit. For the Sommer scale at the physical 
point, a naive linear extrapolation to the continuum limit 
gives r = 0.597(24) (±?|) and 0.612(33)(±^) fm for 7 * /G 
using the PS and V meson dispersion relations, respec- 
tively. The systematic errors are estimated comparing 
the results of various combinations of chiral extrapola- 
tions (linear and quadratic fits for my and ro). We find 
that the results using 7jSy G (PS) and 7^/g(V) are consis- 
tent in the continuum limit within the statistical errors. 
A constrained fit requiring the same continuum value, as 
shown in Fig.|H| leads to r = 0.603(19) (j^) fm, where 
the statistical error was estimated neglecting the corre- 
lation between the PS and V results and the systematic 
error was estimated from the results of constrained fits 
using various combinations of chiral extrapolations for 
my and ro- 



FIG. 20: Sommer scale at the physical point as a function 
of the lattice spacing. Errors are statistical. Lines show the 
continuum extrapolation with a constraint requiring the same 
continuum value for PS and V results. 



to temporal Wilson loops. In order to enhance the over- 
lap C(r) with the ground state, we smear spatial links 
[3lj . We fit the data at t « 0.45 - 0.90 fm where a large 
overlap is observed. As in previous studies on isotropic 
lattices, we find no apparent string breaking effects in 
V(r). Therefore, we fit the potential with 



V(r) = A+ - 
r 



(34) 



for the range r m 0.35 — 1.5 fm. Parameters for the 
potential calculations are summarized in Table IVT1 
With J2U, the Sommer scale ro 1^2 defined by 



,dV(r) 



dr 



= 1.65 



(35) 



D. Beta functions 

Finally we attempt a rough estimation of beta func- 
tions 



8k 
8a s 



8(3 
8a s 



(37) 



along lines of constant physics defined by mpg/my — 
constant. These quantities are required in a calculation of 
the equation of state in thermal QCD [2(j • We calculate 
the beta functions by 



Ok 



d(m v a t ) 



d(mvo.t 



(38) 



d(mps/mv) d(mps/mv) 

( d(myat) d(m PS /my) 
dp 80 

d(myat) d(mps /my) 
dn Ok , 



usin g th e data for m ps and my listed in Table IIIII In 
Ref. [20] , a slightly different method was adopted because 
the matrix in the right hand side of (|38|l sometimes be- 
comes almost singular in the large quark mass region. 
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P 

FIG. 21: Beta functions at 7^ /G (PS) for £ = 2. Thick 
curves are the results for given values of mps/mv, while thin 
curves represent their errors. 



Since quarks are not quite heavy in this study, we adopt 
the simpler method using i|38|) . We fit mps and my to 
the general quadratic ansatz in (3 and k. Because the 
data noticeably deviates from the quadratic form when 
we include all values of k, we restrict ourselves to three 
k's around the target mps/mv in the fit, while all four 
/3's are included. Our estimates for the beta functions 
are summarized in Figs. 1211 and 1221 



VII. CONCLUSIONS 

In this article we initiated a systematic study of two- 
flavor full QCD on anisotropic lattices. We have deter- 
mined, for the clover-improved Wilson quarks coupled 
to an RG-improved glue, the bare anisotropy parame- 
ters which realize a consistent renormalized anisotropy 
£ = 2 in both quark and gauge sectors. In the quark 
sector we employed both pseudo scalar and vector me- 
son channels for calibration. The results for the bare 
anisotropy parameters are summarized in Eqs. (|25[) - l|30|l 
as functions of j3 and k for the range a s « 0.28-0.16 fm 
and mps/my ~ 0.6-0.9. The difference between the two 



-0.02 1 ' ' ' ' 1 ' ' ' ' 1 ' ' ' ' 1 

1.8 1.9 2.0 2.1 

P 

FIG. 22: Same as Fig.|2T]but at 7^ /G (V). 



calibration methods should be 0(a). We confirmed that 
the difference in the bare anisotropy parameters actually 
vanishes toward the continuum limit. 

We have also attempted to calculate some basic quan- 
tities using data measured in the runs made for the cali- 
bration and interpolating them to the point correspond- 
ing to £ = 2. Although errors from interpolations are 
introduced, this enabled us to carry out an initial deter- 
mination of the lattice scale and beta functions. For the 
Sommer scale r , we found that r from different calibra- 
tion methods led to a consistent value in the continuum 
limit. 

We wish to apply our results to study of heavy quarks 
and thermal QCD, in which simulations can be di- 
rectly made with £ = 2 anisotropic lattices using the 
parametrizations Eqs. (|25Jl -130 (l . We hope to report on 
such studies in the near future. 

This work is in part supported by the Large-Scale 
Numerical Simulation Project of the Science Informa- 
tion Processing Center (S.I.P.C) of the University of 
Tsukuba and by the Grants-in-Aid for Scientific Re- 
search by the ministry of Education (Nos. 12304011, 
12640253, 12740133, 13640259, 13640260, 13135204, 
14046202, 11640294). Simulations are carried out on Fu- 



15 



0.2 




0.6 0.7 0.8 0.9 

m ps /m v 



FIG. 23: ?F ass ratio (PS)-^ isp (PS) as a function oimps/mv 
at various values of j3 and re. 

jitsu VPP5000 at S.I.P.C. 



Appendix A: £p from the ratio of screening and 
temporal masses 

In this paper, we adopt the dispersion relation for 
mesons for defining the fermionic anisotropy £p. An 
alternative definition of £p is given by the ratio of the 
masses measures in a spatial direction (screening mass 
m s ) and the temporal direction (temporal mass mt), 

^mass ratio _ / _ /on\ 

= m s /m t , (39) 

as adopted in a quenched study for PS mesons For 
clarity, we denote £p defined by (|18fl using the dispersion 



relation as £p I&p ' in this appendix. 

A disadvantage of the mass ratio method is that, to 
obtain reliable values of m s and m t suppressing con- 
tamination of excited states, we need to prepare well- 
tuned smeared sources both in the spatial and tempo- 
ral directions, and/or carry out multi-pole fits, on suffi- 
ciently large lattices. In this paper, because we do not 
have propagators with temporally smeared sources, we 
study spatial propagators with point-point source and 
sink. We find that, when quarks are light, the effec- 
tive mass of the spatial PS meson correlator does not 
show a clear plateau, and sometimes shows a decreasing 
tendency even at the maximum distance x = N s /2 — 1. 
This means that our spatial lattice sizes 8-12 may not be 
large enough to suppress excited states. Unfortunately, 
the number of data points is also not sufficient to at- 
tempt a multi-pole fit. Therefore, in the following, we 
just adopt the value of effective mass at x = N s /2 — 1 for- 
ms. Strictly speaking, this value gives an upper bound 
of m s . Therefore, the resulting £™ ass ratl ° may be larger 
than the true value when quarks are light. 

Because our temporal lattice sizes are sufficiently large, 
we do not encounter a similar problem in the calculations 
of m t and t; F lsp ' ■ 

Figure I2H1 shows the difference of the two fermionic 
anisotropics ££? ass ratio (PS) - ^ isp (PS) as a function of 
mps/m v at j3 = 1.9, 2.0 and 2.1. We find that the dif- 
ference is consistent with zero when mps/my is small. 
A slight overshooting at mps/ray <0.7 may be under- 
stood by the fact that our £™ ass ratl ° is an upper bound 
for the true value as discussed above. We also find that 
the difference at large quark masses decreases toward the 
continuum limit. At (3 >2.0 (2.1), two methods are con- 
sistent with each other at nips/my ^0.75 (0.8). 
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